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Diffraction phenomena of large molecules have been studied in many experiments, and these 
experiments are described by many theoretical works. In this paper, we study Ceo molecules 
single and double-slit diffraction with quantum theory approach, and we pay close attention to 
the Ceo diffraction experiment carried out by A.Zeilinger et.at in 1999. In double-slit diffraction, 
we consider the decoherence effect, and find the theoretical results are good agreement with 
experimental data. 
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1. Introduction 

As is well known, the matter-wave diffraction has become a large field of interest over the last years, 
and it is extended to electron, neutron, Cgo, atom, more massive, complex objects, like large molecules 
-^2, C'eo a-nd C70, which were found in experiments [1-5]. At present. There are classical and quantum 
methods to study interference and diffraction [6-12]. The classical optics with its standard wave-theoretical 
methods and approximations, in particular those of Huygens and Kirchhoff, has been successfully applied to 
classical optics, and has yielded good agreement with many experiments. This simple wave-optical approach 
also gives a description of matter wave diffraction. However, matter-wave interference and diffraction are 
quantum phenomena, and its full description needs quantum mechanical approach. Recently, there are some 
quantum theory approach to study electron and Cgo diffraction, and obtain some important and new results 
[13-17]. In this paper, we study the Cgo single and double-slit diffraction with the quantum approach, and 
compare to the Cgo diffraction (A « 250A) carried out by A. Zeilinger et. at in 1999 [18]. The results we 
obtained are shown in Fig. 3, Fig. 4 and Fig. 5. In view of quantum mechanics, the Cgo has the nature 
of wave, and the wave is described by wave function il}{r,t)^ and the wave function ip{r,t) has statistical 
meaning, i.e., | iplj^t) ^ can be explained as particle's probability density. For the Cgo slit diffraction, if 
we can calculate the Ceo wave function i^{f,t) distributing on display screen, we can obtain the diffraction 
intensity of Ceo molecules, since the diffraction intensity is directly proportional to | ip{r,t) p. In the slit 
diffraction, the Ceo wave functions can be divided into three parts. The first is the incident area, and the Ceo 
wave function is a plane wave. The second is the slit area, where the Ceo wave function can be calculated 
by the Schrodinger wave equation. The third is the diffraction area, where the Ceo wave function can be 
obtained by path integral. In the following, we shall calculate these wave functions. Finally, we study 
the decoherence effect in the Ceo double slit diffraction. The decoherence effect is important for the large 
molecules diffraction, and the result is agreement with the experiment data. 
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FIG. 1: Double-slit geometry with a the slit width, b the slit length and d the distance between the two slits. 



2. Quantum approach of Ceo diffraction 

In an infinite plane, we consider a double-slit, its width a, length b and the slit-to-slit distance d are shown 
in FIG. 1. The x axis is along the slit length b and the y axis is along the slit width a. We calculate the Cgo 
wave function in the first single slit (left) with the Schrodingcr equation, and the Cqo wave hmction of the 
second single slit (right) can be obtained easily. At time t, we suppose that the incident plane wave travels 
along the z axis. It is 

Mz,t) = Aei^P^-^'\ (1) 

where A is plane wave amplitude. 

The potential in the first single slit is 

V(xvz) = (^ 0<^<^0<y<a,0<^<c, 
^ ' ^' ' 1^ 00 otherwise, ^ ' 

where c is the thickness of the single slit. The time-dependent and time-independent Schrodinger equations 
are 

d fl^ Q2 

^ + ^ + ^-^«^1=0. (4) 

where M{E) is the mass(energy) of the Ceo. The relation between ilj{r,t) and is 

il>{fi,f2---rN,t) =ip{fi,f2---fN)9{t). (5) 

The Eq. (19) become as 

In Eq. (4), the wave function ip{x, y, z) satisfies the boundary conditions 

i;{0,y,z) = i;{b,y,z) = 0, (6) 

tp{x,0,z) = 'ip{x,ai,z) = 0. (7) 
The partial differential Eq. (4) can be solved by the method of separation of variable. By writing 



i;{x,y,z) = X{x)Y{y)Z{z). 



(8) 



The general solution of Eq. (3) is 

tlJi{x,y,z,t) = '^ipmn{x,y,z,t) 
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Equation (9) is the Cgo wave function in the first single slit. Since the wave functions are continuous at 
z = 0, we have 

il>o{x,y,z,t) \z=o= ipiix, y,z,t) \z=o, (10) 
from Eqs. (2), (6) and (9), we can obtain the Fourier coefficient -Dmn by Fourier transform 
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otherwise, 
substituting Eq. (11) into Eq. (9), we can obtain the Ceo wave function in the first single slit. 
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(11) 
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The Ceo wave function in the second single slit can be obtained by making the coordinate translations 
x' = X, y' = y — a — d, z' = z, and we can obtain the Ceo wave function ^2(2;, y, z, t) in the second slit 

, 16A . (2n+l)7rx . (2m + l)7r(y - ai - d) 

■4>2{x,y,z,t) = > — — sm^ — sm '-^ 

^ 2m + 1 (2n + llTT^ a 

.e*V-p- — ^ "^e-ft^*. (13) 

3. The wave function of Ceo diffraction 

With the approach of path integral, we can calculate Ceo wave function in the diffraction area. 

The diffraction area is shown in FIG. 2, where ro is the position of point P(j{xo,yoT c) on the slit surface 
(z = c), P{x,y,z) is an arbitrary point in the diffraction area. From Eq. (12), we can obtain the wave 
function in slit surface {z = c) 

, , ^ 16^ . (2n + l)7ra;o . (2m+l)7r2/o 

ipP^{x,y,z,t) = > — — — — — — ^sm sm 

(2m + l)(2n + 1 TT^ a 

m,n— 



• e'V-p- — p '=e""R^*°. (14) 

The diffraction wave function ^out on the point P{x, y, z) can be calculated by the formula of path integral 
[17] 

tpp{r,t) = j K{r,t;ro,to)tpPa{ro,to)dxodyo, (15) 
where K{r,t;ro,to) is the Ceo molecules propagater, it is 
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FIG. 2: Diffraction area of the single slit 



where M is the mass of Ceo molecules, and R is the distance between point Pq and point P, it is 



R = ^/\r^^r^ =^ V(x - xo)2 + (y - ?/o)^ + - c)2 

= + 2/2 _|_ ^2 _|_ 3,2 _|_ y2 _|_ ^2 _ 2a;2;(, - 2yyo - 2zc. (17) 

Assume that the angle between r and x axis (y axis) is J — a (J — /3), and then a (/3) is the angle between 
r and the surface of yz{xz). we obtain 

R^ = + + xl + yl + - 2xxo - 2yyo - 2zc 

= — 2r sin axo — 2r sin pyo + x^+y^ 

~ _ 2rsinaa;o — 2rsmpyo, (18) 

where = x"^ + y'^ + {z — c)^, In Eq. (18), the terms Xq and j/q can be neglected. Substituting Eqs. (14), 
(16), (18) into (15) yields 
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With de Broglie relationship p = hk = y, there is 
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Equation (24) is tlie Cgo diffraction wave function of the first slit, and the diffraction wave function ip2P 
for the second slit can be obtained by making the coordinate translations x' = x,y' = y — {a + d), z' — z, it 
is 

m— n— ^ 



2ME (2n + l)7r (2m+l)^ 2 ^ " ^2l 

- ( r - ( : + i^k - -5) Wcos2 a - (-)2] 



I 



'e-^^«'"«-'sini^I^x'rfx' 
b 

1a+d 



f " g-ifesin/3V sin _ ^a + d))dy', (25) 

where d is the two slit distance. The total diffraction wave function for the double-slit is 

V'p(a;, y, z, t) = ciipip{x, y, z, t) + C2il>2p{x, y, z, t), (26) 

where ci and C2 are superposition coefficients , and |cip + |c2|^ = 1. For the Ceo single-slit diffraction, we 
can obtain the relative diffraction intensity / on the display screen. 



I ^\Mx,y,z,t)\''. (27) 
For the Ceo double-slit diffraction, we can obtain the relative diffraction intensity I on the display screen, 
I oc \il)p{x,y,z,t)f' 

= cf |V'ip(a;, y, z, t)\^ + c||^2p(a;, y, z, t)\^ + 2ciC2i?e[V'tp(a;, y, z, i)V'2p(a;, y, z, t)]. (28) 

From Eqs. (27) and (28), we can obtain the relation between diffraction intensity and diffraction angle. In 
Ref. [16], their experiment data is about the relation between diffraction intensity and diffraction position. 
The relation is 

sin[l3\ = ^= ^ (29) 

Where I is the distance from the slit to display screen, and s is diffraction position. 
4. Decoherence effect in double-slit diffraction 

Decoherence is introduced here using a simple phenomenological theoretical model that assumes an ex- 
ponential damping of the interferences [7, 20, 21], i.e., the decoherence is the dynamic suppression of the 
interference terms owing to the interaction between system and environment. The Eq. (26) describes the 
coherence state coherence superposition, without considering the interaction of system with external en- 
vironment. When we consider the effect of external environment, the total wave function of system and 
environment for the double-slit factorizes as [7] 

■tpoutix, y, z, t) = ciipouti O \Ei >t +C2ilJout2 O \E2 >t, (30) 

where \Ei >t and \E2 >t describe the state of the environment. The diffraction intensity on the screen is 
now given by [7]: 

I = {1 + \at\'^){cl\'ipoun{x,y,z,t)\'^ + cl\jpout2{x,y,z,t)\'^ + 2ciC2AtRe[ijl^f^{x,y,z^ (31) 

where at =t< -E2I-B1 >t, and = . Thus, At is defined as the quantum coherence degree. In Eq. 

(31), the two slits wave functions ipip and ip2P are calculated by the quantum approach (in Eqs. (24)-(25)). 
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In Rcfs. [7], the two slits wave functions are two Gaussian wave packets. The fringe visibility of v is defined 
as [7]: 

(32) 



+ In 

Where Imax and are the intensities corresponding to the central maximum and the first minimum next 
to it, respectively. The value for the fringe visibility oi v = 0.50 is obtained in Zeilinger et. al. experiment 
[18] {Imax = 880, /min = 300), and the quantum coherence degree At « [7]. 

5. Numerical result 

Next, we present our nimierical calculation of relative diffraction intensity. The main input parameters 
are: Ceo molecules mass M = 1.4 x 10~^^kg, and Planck's constant h = 1.055 x 10~^^Js. For the single-slit 
experiment [18], the Ceo velocity v = 220m/ s (corresponding to Ceo wave length A = 250A), the slit width 
a = 10/xm, and the distance between slit and display screen / = 2.29m. In our calculation, we take the 
same experiment parameters above, and the theoretical input amplitude parameters are: A = 2.87 x 10^^, 
the diffraction angle on yz surface a = rad, the slit length b = 0.01m and the slit thickness c = 1.3/zm. 
From Eq. (27), wc can obtain the single slit diffraction intensity pattern, and it is shown in FIG. 3. In 
FIG. 3, the solid curve is our calculation result, and the dot curve is the experiment data [18]. From FIG. 
3, we can find the calculation result is agreement with experiment data. For the double-slit diffraction, we 
consider two cases: coherence superposition and decoherence effect. For coherence superposition, from Eq. 
(28), we can calculate the diffraction intensity and it is shown in FIG. 4. The experiment parameters are: 
the Ceo wavelength A = 250A, the first and second slit width, a = O.Obfim, the distance between the two 
slit d = 0.05/im, the slit thickness c = 1.3/um, the distance between slit and display screen / = 1.25m, In 
our calculation, we take the same experiment parameters above, and the theoretical input parameters are 
superposition coefficients ci = 0.566, C2 = 0.824 (|cip + |c2p = 1) and amplitude parameter A = 1.27 x 10^^. 
In FIG. 4, the solid curve is our theoretical calculation, and the dot curve is the experiment data [18]. From 
the FIG. 4, we find that the theoretical result is in accordance with the experiment data, when the position 
s is in the range of \s\ > 100/um. When the position s is in the range of \s\ < 100/xm, the theoretical result 
has a large discrepancy with the experiment data. We find the discrepancy can be eliminated when the 
decoherence effect is considered. Prom Eq. (31), we can obtain the diffraction intensity pattern and it is 
shown in PIG. 5. In calculation, superposition coefficients ci = 0.565, C2 = 0.824, amplitude A= 1.69 x 10^^ 
and quantum coherence degree At = 0.50. From PIG. 5, we can find that the new calculation result is more 
improvement than it in FIG. 4, and it is in accordance with the experiment data, i.e., when the decoherence 
effect is considered, the discrepancy between the theoretical result and experiment data can be eliminated. 

6. Conclusion 

In conclusion, we study Ceo single and double-slit diffraction with quantum theory approach. The 
theoretical result of single-slit is in accordance with the experiment data. For the double-slit diffraction, we 
study the diffraction intensity by two approaches, which are the coherence superposition and decoherence 
mechanism. When we only consider the superposition of coherence, we find the theoretical result has a 
large discrepancy with the experiment data, and when we consider the decoherence mechanism, we find the 
theoretical result is in accordance with the experiment data. So, the decoherence mechanism is important 
for the double-slit diffraction of the large molecules. We think the new quantum theory approach has 
universal applicability, such as, it can study electron, atom and molecular diffraction diffraction, and it can 
also be studied multi-slit and grating diffraction. 
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FIG. 3: Comparison between theoretical prediction from Eq. (27) (solid line) and 
experimental data taken from [18] (circle point) for Ceo single-slit diffraction. 
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FIG. 4: Comparison between theoretical prediction from Eq. (28) (solid line) and experimental data taken 
from [18] (circle point) for Ceo double-slit diffraction, no including the decoherence effects 




FIG. 5: Comparison between theoretical prediction from Eq. (31) (solid line) and experimental data taken 
from [18] (circle point) for Ceo double-slit diffraction, including the decoherence effects. 



